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cos6-71-4 The study of the anisotropy of the arrival directions is an essential tool to investigate 
the origin and propagation of cosmic rays primaries. A simple way of recording many cosmic rays 
is to record coincidences between a number of detectors. We have monitored multi-TeV cosmic rays 
by a small array of water cherenkov detectors in Tehran(35°43' N, 51°20'E, 1200m a.s.l). More than 
1.1 X 10® extensive air shower events were recorded. In addition to the Compton- Getting effect 
due to the motion of the earth in the Galaxy, an anisotropy has been observed which is due to a 
' unidirectional anisotropy of cosmic ray flow along the Galactic arms. 
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I. INTRODUCTION 



Although Cosmic Rays( CRs) have been known for almost one century, their origin remains uncertain, mostly 
0^ ■ because their trajectories are bent by Galactic magnetic fields and they do not individually point back to their 
C^l \ sources. Cosmic rays in the lower energy range have gyro radii of about Ipc or less in typical galactic magnetic fields( 
a proton with an energy of lO^^eV would have a gyro radius of 1 pc in a 1/iG field). Moreover, since these fields are 
' (— I '' chaotic on scales ranging at least from lO^cm to lO^^cm ll!], the transport of CRs is diffusive up to high energies, which 
^ i: tends to make their angular distribution isotropic. Therefore, even collectively, the CR arrival directions hold virtually 
no information about the source distribution in space. However, as the energy of the CRs increases, it acan appear 
either because the diffusive approximation does not hold anymore, or because the diffusion coefficient becomes large 
enough to reveal intrinsic inhomogeneities in the source distribution. Specifically, even if the diffusive regime holds, 
C3 '. the density of CR sources in the Galaxy is believed to be larger in the inner regions than in the outer ones, and this 
can cause a slightly higher CR flux coming from the Galactic center than from the anti-center. Meanwhile, the global 
I [ CR streaming away from the Galactic plane( towards the halo) can be a source of measurable anisotropy. However, 
^ . the detailed angular distribution of CRs is quite hard to predict, even if we assume a definite source distribution, 
' because it also depends on the propagation conditions, which are related to both large scale and small scale magnetic 
. field configurations, and on the position of the Earth relative to major magnetic structures, such as the local Galactic 
C^l ' arm. 

From a general point of view, the characterization of the CR anisotropy provides useful information to constrain the 
GCR diffusion models, notably the effective diffusion coefficients, related to the magnetic field structure. Indeed the 
\ level of CR anisotropy depends on the diffusion coefficient, £>: in a simple model where CR sources are homogeneously 
C . distributed in a disk of thickness 2h and the CRs are confined in a halo of height H, the anisotropy at a distance z 
■ above the Galactic plane(z < h) is estimated as A ~ 3D/cH x z/h The deviation from isotropy is typically below 
^ , 1% and can be as low as 0.03% Anisotropy measurements at various energies can thus provide crucial information 
' about the energy dependence of the diffusion coefficient. This information is particularly important to constrain the 
r\ I GCR source spectrum, since it sets the relation between the source power law index and the observed one, through 
• the energy dependent confinement of CRs in the Galaxy. This diffusion maybe is broadly along the magnetic field 
lines which are in tubes of dimensions greater than the gyro radii. So the direction of the peak of the anisotropy 
would indicate the direction back towards the cosmic ray source, and the amplitude of the anisotropy would give 
information on the scattering process involved in the diffusion. Specifically, an estimate of the mean free path might 
be obtained. 

The anisotropy is due to a combination of effects. Compton and Getting^] proposed in 1935 that the motion of the 
solar system relative to the rest frame of the cosmic ray plasma should cause an energy independent dipole anisotropy 
with maximum in the direction of motion. The earth's rapid motion in space, resulting from the rotation of our 
galaxy, results variations in cosmic ray intensity fore and aft of the earth's motion. Following Compton - Getting 
(1935), the magnitude of the anisotropy is expressed as 

?7= (7 + 2)-cos/?, (1) 
c 
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where 7 denotes the power law index of the energy spectrum of cosmic rays, u the velocity of the detector relative 
to the production frame of the cosmic rays ( where they are presumed to be isotropic), c the speed of light, and /? 
the cosmic ray direction relative to u, i.e. cos/3 is the projection of the cosmic ray along the forward direction of u. 
In fact the value of (7 + 2)^ is {fmax - frmn)/ifmax + frnin) with fmax thc countiug rate along the direction of the 
velocity and fmin along the contrary direction. The magnitude of the anisotropy is extremely small and independent 
of the cosmic ray energy. Our data will be analyzed in a sun-centered frame, and so if data accumulation is done 
over an integer number of solar years, it is only necessary that the orbital speed of the Earth around thc sun(~ 30km 
s"""^) be considered. The large effect due to the Galactic rotational speed (220 km s~^) will cancel out as the data are 
averaged over this time @. Many experiments have been carried out for detection of this effect 0, Q . 

Doppler effect studies of globular clusters and extra galactic nebulae have revealed a motion of the earth of about 
220km s""'^ towards right ascension a ~ 21h and declination S ~ 47°N due chiefly to the rotation of the Galaxy. 
This motion, with a speed of about 0.1% c will affect the intensity of the incoming cosmic rays by changing both 
the energy of the cosmic ray particles and the number received per second. Using value of 220 kms^^ for u, and 2.7 
for the spectral index, Eqn. (1) gives a Compton Getting Effect(CGE) amplitude of 0.345 x 10^^ for the fractional 
forward-backward asymmetry caused by the motion of thc earth in thc Galaxy. 

Another effect that can produce sidereal modulation is solar diurnal and seasonal changes in the atmospheric 
temperature and pressure. As the atmospheric temperature and pressure change during the course of a day, the 
balance of cosmic ray secondary particle interaction and decay changes. This propagates to changes in the detection 
rate that depend on the detector type( air shower, underground muon, surface muon) and on the energy threshold. 
These changes tend to have a strong Fourier component with a frequency of one solar day(~ 365 cycles/ year) and one 
solar year( 1 cycle/ year). In some( but by no means all) experiments, the interplay between the daily and seasonal 
modulation can produce significant modulation in sideband frequencies of ~ (365± l)cycles/ yearQ. The modulation 
with frequency 366 cycles/ year appears as a sidereal modulation. The size of the atmospheric contribution to apparent 
sidereal anisotropy can be estimated from the amplitude of the pseudo- sidereal( 365 cycles/ year) frequency. If it 
is large, the atmospheric effect can be subtracted using the amplitude and phase of the pseudo- sidereal component. 
The anisotropy that remains after accounting for the Compton- Getting and atmospheric effects is due to solar and 
galactic effects. At the lowest energies(~ lOOGeV), the Interplanetary Magnetic Field( IMF) produced by the solar 
wind effects the sidereal anisotropy: when the local IMF points toward the sun, thc anisotropy peaks at about 18hr 
right ascension, while it peaks at 6hr when it points awayfs', The average over the two configurations produces 
a small, residual anisotropy peaking at around 2-4hr. At higher energies, local IMF plays a negligible role. Instead, 
the heliosphere extending to distances of order lOOAU is believed to induce anisotropy in cosmic rays with energies 
around lTeVp^.[TH. Beyond this energy, the anisotropy is believed to be primarily of galactic origin. For instance, the 
galactic magnetic field around the solar system neighborhood could produce anisotropy. Also, an uneven distribution 
of sources of cosmic rays ( presumably, mostly supernova remnants) may produce anisotropy. It is believed that star 
formation( and thus, supernova remnants) occur primarily in the spiral arms of the galaxy. The earth is located at 
the inner edge of the Orion spur. Thus, in the direction of the Orion spur{ galactic longitude between 60° to 270°) 
they are distributed nearby sources of cosmic rays, while in the complementary direction, they are much further away. 

Because of small anisotropy, large data sets are required to make useful measurements which overcome the statis- 
tical uncertainties of counting experiments. A simple way of recording many cosmic rays is to record coincidences 
between a number of detectors. Few statistically significant anisotropics have been reported from extensive air shower 
experiments in the two last decades. Aglietta et al.( 1996,EAS-T0P)[12J published an amplitude of (3.7 ±0.6) x 10"'* 
and phase — (1.8 ± 0.5)hr LST, at Eo ~ 200TeV. Analysing the Akeno experiment, Kifune et al.(1986)^| re- 
ported results of about 2 x 10~^ at about 5 to 10 PeV. An overview of experimental results can be found in |14|.We 
have operated a small array of water cherenkov detectors on the roof of physics department at Sharif University of 
Technology in Tehran(35°43' N, 51°20'E, 1200m a. s. 1=890 gcm~^) as a prototype for constructing an Extensive Air 
Shower(EAS) array on Alborz mountain range at altitude of over 2500m near Tehran. 

The main purpose of this article is to study the unidirectional anisotropy of cosmic rays flow along the Galactic 
arms which was observed in the sidereal time at energies above 50 TeV. We describe the experimental setup in section 
II, and the data analysis and discussion in section III. 

II. EXPERIMENTAL SETUP 

Four Water Cherenkov Detectors( WCDs) used for recording Extensive Air Showers( EASs). The four WCDs 
consist of cylindrical tanks made of polyethylene with a diameter of 64cm and a height of 130cm filled up to a height 
of 120cm with 382 liters of purified water. All the inner surfaces of the four Cherenkov tanks were optically sealed and 
covered with white paint which reflects light in a diffusive way. Each one of them have a single 5.2cm PMT( model 
EMI 9813KB) located at the top of the water level along the cylinder axis. The array arranged in a square with side 
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608cm as shown in Fig.l, on the roof of the Physics Department at Sharif university of Technology in Tehran(35°43' 
N, 51°20'E, 1200m a. s. 1=890 gcm~^). The signal produced by secondary particles of an EAS are triggered with an 
amplitude threshold of -500mV by an 8-fold fast discriminator ( CAEN N413A). The threshold of each discriminator is 
set at the separation point between the signal and background noise levels. The discriminator outputs are connected 
to one of three Time to Amplitude Converters( TAG) (EG&G ORTEC 566) which are set to a fuU scale of 200ns ( 
maximum acceptable time difference between two WGDs). The output of due to the No. 3 WCD is connected to the 
start inputs of TAGl, TAC2, and TAC3. The outputs of due to No.l, 2, and 4 WCDs are respectively connected to the 
stop inputs of TAGl, TAG2, and TAC3. then the outputs of these three TAGs are fed into a multi-parameter Multi 
channel Analyzer( MCA) ( KIAN AFROUZ Inc.) via an Analogue to Digital Converter ( ADC)( KIAN AFROUZ 
Inc.) unit. The output of TACl triggers the ADC, and the 3 time lags between the output signals of PMTs (3,1), 
(3,2), and (3,4) are read out as parameters 1 to 3. So by this procedure an event is logged. 

III. DATA ANALYSIS AND DISCUSSION 
A. Array event rate 

The data set covers a total period of 2.7 x 10^ seconds. A total of 1.1 x 10^ events with zenith angles < 60° 
were collected during this time, giving a mean event rate of one event every 24.5 seconds. Figure 2 shows the events 
time- spacing distribution. Since events arrive randomly in time, it is expected that this will follow an exponential 
distribution, viz. 

/(0-/(0)cxp(-i/T). (2) 

The event rate can be obtained by fitting this function on the events time-spacing distribution. One event per every 
T = 24.1s is obtained from the fit. A non-random component for the cosmic ray flux for example a point source of 
gamma ray gives rise to deviation from the exponential law. Our observed distribution is in good agreement with the 
exponential law. 

B. Atmospheric effects on counting rate 

The rate of shower detections depends on a number of factors. If either temperature or pressure variations have 
Fourier components in solar or sidereal times, spurious components may be introduced into the shower detection 
rate 0] . Various methods are used in order to study the dependence of event rate on atmospheric ground pressure, 
p, and temperature, T, p^ . In a multiple regression analysis of event rate against pressure and temperature, the 
temperature effect is found to be statistically insignificant. The CR intensity dependence on barometric pressures in 
half- hour intervals are shown in Fig. 3. We can describe the dependence by the following function 

i? = i?ocxp( ^°"^' ) (3) 
Pi 

The values of Rg = 74 events per half an hour, po=883.7mb, and pi = 171.5mb were obtained from the data. Pi denotes 
the measured air pressure at a given time. By this empirical function, we weighted the raw event rates for atmospheric 
pressure. Figure 4 shows the mean half-hour event rate distributions with and without correction for atmospheric 
ground pressure. The distribution of corrected rates is consistent with a Gaussian distribution, as expected for the 
statistical fluctuation of the events rate, and there is no residual temperature effect. 

C. Zenithal angle distribution of the EAS events 

Since the thickness of the atmosphere increases with increasing zenith angle, 9, the number of EAS events is strongly 
related to 9, as shown in Fig. 5. The differential zenith angle distribution can be represented by 

dN = constant. {5i cos 6 + 82 sin 9) cos" 9 sin 9d9. (4) 

Where we split into particles entering through the lid of cylinderical tank of WCD or through its walls. The first 
term in parenthesis of Eqn.(4) is related to the lid and the second to the walls. The parameter 5i includes the area 
of lid surface, S'l, and detection probability of particles entering through the lid. Pi. The parameter 82 also includes 
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the greatest surface area of vertical profile of the WCD, S2, and detection probability of particles entering through 
its walls, P2. So we can split Sj{j — 1, 2) in the form of 



Sj^S.P,, (5) 

where only Pjs are determined from the simulation ( see the Appendix). 5*1 and S2 are respectively 3.2 x lO'^cm^ and 
7.68 X lO^cm^. By fitting Eqn (4) on our experimental data, n — 7.3 is obtained. 



D. Energy Threshold of our experiment 



Since we can not determine the energy of the showers on an event-by-event basis, we estimate the energy threshold 
of our array by the CORSIKA code for simulation of EAS events 16]. In order to record a shower it is necessary 
at least one particle passes through each of the four WCDs. Because our array has been arranged in a square with 
side 608cm we can detect a shower if density of secondary particles to be at least pr = 1/Seff particles cm~^ at 
r — 608 X \/2/2 = 420cm, where Seff is the effective surface area of a WCD. This area is calculated as follows: 



J^^^iPiSi cose + P2S2 sine) smedO 



Seff = ^ .,,3 = 6.5 X 10^ cm^ (6) 



with Pi = 0.88, Si = 3.2 x lO^cm^, P2 = 0.93, and S2 = 7.68 x lO^cm^. The upper limit, 7r/3, is due to events 
selection with zenith angles < 7r/3. We simulated more than 10^ EAS events with CORSIKA code using the hadronic 
interaction models QGSJET and GHEISHA. The energy range for primary particles was selected from 5TeV to 5PeV, 
with differential flux given by dN/dE (x E^'^^'' . These simulations are in different directions with zenith angle from 
0° to 60° and azimuth angle from 0° to 360°. Finally from the CORSIKA simulations with — 1.5 x 10~'*particles 
cm~^ at r = 420cm, we obtained the energy threshold Eth = 50TeV. 



E. Sidereal time distribution 



After atmospheric correction we calculated the sidereal time (ST) from ST = STo+a{ZT— ZTq) . STo can be looked 
up in an almanac fv^ for the time ZTq, ZT is the solar time, and a — 1.00273790935. Figure 6 shows percentage 
variation in intensity of the cosmic rays with sidereal time. The data have been fit to Eqn. (7) which describes a curve 
with first and second harmonics( i.e with a once-per-day and a twice-per-day variation), 

y = Ai cos[^(t - Ti)] + An cos[^(t - Tn)]. (7) 

Where t is in hours. The fitting results of data are summarized in Table 1. The CGE would contribute to the component 
Aj in the sidereal time asymmetry. This analysis shows that the anisotropy has a peak close to the sidereal time 21h, 
when the zenith is toward the earth's motion. The amplitude of the first harmonic is 0.32%. So there is a definite 
sidereal time variation whose phase and amplitude are close to those predicted. In order to calculate the magnitude 
of anisotropy due to CGE i.e the value rj in Eqn.(l), a mean value for cos/3 is needed. Assume S is the declination of 
the direction of earth's motion, A the latitude of the observer and H the hour angle between the observer's meridian 
and the direction of motion, then to consider Fig. 7 the angle <j> between the observer's zenith and the direction of 
earth's motion is given by 

cos = sin 5 sin A -f cos (5 cos A cos i?. (8) 

On the other hand, cos (3 is calculated by 

cos/3 = cos(/)cos6' -I- sin^sinScoso;. (9) 

Where 9 is the zenith angle of cosmic ray and a difference between the azimuth angle of the direction of motion and 
of cosmic ray( Fig. 7), that is a = Ai — A2, which Ai and A2 are obtained by 



sin S — sin A cos cj) + cos A sin (j) cos Ai , 



(10) 
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and 

sin (5' = sin Acos6' + cos A sin 6* cos A2 . (11) 

Where 6' is the dechnation of cosmic ray. According to Eqns (8) to (11), the 24-hour mean of the component of cosmic 
ray in the direction of motion (cos/3) may be obtained. Using A = 35°43 and 6 = 47°, we calculated the 24-hour 
mean value of coscp — 0.43 with Eqn.(8). With distribution of {PiSi cos9 + P2S2 sm9) cos^ "^ ^sin^ which describes 
the acceptance of detectors and the cosmic ray absorption in the Earth's atmosphere because of inclination from the 
vertical direction in Tehran ( section III, C) we calculated the mean value of cos6'=:0.88. Also the mean values of Ai 
and A2 were obtained by using the mean value 5' . Figure 8 shows the distribution of cosmic rays declination. The 
mean value of declination is 5' = 32.5°. From Eqns. (10) and (11), Ai and A2 obtain 49° and 86° respectively. Finally 
from Eqn.(9) a value of 0.72 is obtained for cos/3 and this is multiplied by the expected CGE amplitude of 0.345% 
to yield a predicted effect of expected value of 0.248%. The value obtained from experimental data is 0.32% which 
about 0.07% is more than the CGE value. This remaining asymmetry of 0.07%, presumably has an origin different 
to the that of the CGE. 

Since the recorded data are in Tehran with latitude 35°43'A^, the majority of cosmic rays are from the spiral arm 
inwards direction, which is at about 20h in right ascension and 35° in declination [18]. So the remaining asymmetry is 
probably due to unidirectional anisotropy of cosmic ray flow along the Galactic arms. A simple diffusion model [l9| 
suggests that value of this asymmetry, 0.07%, would be roughly equal to the ratio of the scattering mean free path to 
a characteristic dimension of the containment region ( i.e the central Galactic region, with a scale of lOkpc). So with 
amplitude of the anisotropy of 0.07% found in this work, we obtain a mean free path of about 7pc which is about 
perhaps 7 gyro radii. 

Since the anisotropics are low, it is necessary to consider the effect of counting statistics for a finite measured data 
set. If we have N events then the probability of getting a fractional amplitude greater than r is given by [20| . 

P{>r)^e~''°,ko^r^N/A. (12) 

So a convenient parameter for characterizing the anisotropy amplitude probability distribution is /cq- We can take 
V^rrms which corresponds to ko =1, as noise amplitude. For the number of events that we have accumulated, 1.1 x 10^, 
the total amplitude of 0.32% obtained in this work can be arisen by chance with a probability of 0.06 corresponding 
to ko = 2.8. This shows a significant anisotropy ( fco > 1) at the sidereal period. So we conclude that this data set 
gives evidence of anisotropy. 

IV. CONCLUSION 

Cosmic ray data in Alborz observatory clearly shows an anisotropy in sidereal time with the energy threshold of 
^ 50TeV and the mean energy of 121TeV. One part of this anisotropy is due to Earth's motion around the Galaxy 
(the CGE), but our measured asymmetry suggests the possible existence of some other additional effects, probably 
a unidirectional anisotropy of cosmic ray flow along the Galactic arms. The first harmonic amplitude of our total 
measured anisotropy is about 0.32%. The CGE contribution to this anisotropy is about 0.248% and the rest, 0.07%, 
is predicted to be due to the flow along the Galactic arm. The latter anisotropy suggests a mean free path of about 
7pc for these high-energy cosmic rays. The evidence of these anisotropics is based on the value of the parameter feg, 
as suggested by Linsley( 1975) and found in this work to be 2.8, that is, more than fco = 1, the value for the noise 
amplitude. 

The EAS-TOP experiment reported somewhat lower limits in the energy range below 1200TeV 21]. The relatively 
large amplitudes published by the Akeno experiment and our experiment are difficult to reconcile with the results 
of the EAS-TOP experiment. 

Appendix 

In order to obtain the PjS, Eqn.(5), we first calculated the track length of particle passing through lid and walls of 
WCD. To determine the track length distribution we have made the following assumption: 

A) The zenithal and azimuth angles of particles are uniformly distributed. 

B) The random distribution increases linearly with r, the distance from center of lid, i.e. it is proportional to the 
annulus surface 27rrAr. 

The geometry can be simply solved. It is split into particles entering through the lid or through the walls as seen 
in Fig. 1. In these figures, r is the distance from the center of the cylinder, the tank radius is Rq = 32cm, the tank 
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height Hq = 120cm, and ip and arc azimuth and zenith angles of particle, respectively. The simulation process starts 
by randomly choosing r which increases linearly with r. Then ip and are randomly chosen and the track lengths 
within the tank are evaluated by a simple calculation ( the particle could either leave through the wall or the bottom 
lid, as shown in Fig. 1). 

Then the number of photons produced along a flight path was estimated. Charged particles emit light under a 
characteristic angle when passing through a medium if their velocity exceeds the speed of light in the medium. The 
Cherenkov angle is related to the particle velocity and the refractive index of the medium n{ cos0 = l/n(3). For 
relativistic particles, (3=1, and the refractive index of purified water, n = 1.33( for short wavelengths of visible 
region), the Cherenkov light is emitted under 41°. The number of photons produced along a flight path dx in a wave 
length bin dX for a particle carrying unit charge is : 



d^ ~ A2 ' ^ 

where a = 1/137. At wavelengths of 310-470nm the efficiency of the photomultiplier is maximal. Within 1cm flight 
path 220 photons are emitted in this wavelength bin. To consider the effective area of the photomultiplier ( 21.2cm^), 
and neglecting absorption and scattering effects in water we obtained the number of photons received by PMT ( Fig. 
9). Finally with a, rj = 25% quantum efliciency and a G=10** gain for the PMT, the number of electrons prochicod 
in PMT(A^e = NphotonVG) was calculated. As we know, the output signal at the PMT's anode is a current or 
charge pulse. Now to consider the amplitude threshold of discriminator ( -500mV), and the anode load resistor and 
capacitance, we obtained that for producing a pulse with amplitude -500mV, the number of photons received to PMT 
should be more than 1. Because the quantum efficiency of PMT is 25% using Fig. 9 we calculated the detection 
probability, PjS, are respectively 0.88 and 0.93 for lid and walls WCD. 



Acknowledgments 



This research has been partly supported by Grant No. NRCI 1853 of National Research Council of Islamic Republic 
of Iran. 



J.W. Armstrong, B.J. Rickett, and S.R. Spangler, ApJ 443, 209 (1995) 

V. S. Ptuskin, Adv. Space Res. 19, 697 (1997) 

A.G.K. Smith and R.W. Clay, Aust. J. Phys. 50, 827 (1997) 

A.H. Compton, and LA. Getting, Phys. Rev. 47, 817 (1935) 

J. Poirier, and C. D'Andrea, M. Dunford, Proc. 27th ICRC, Hamburg, Germany , 3930 (2001) 

T.F. Lin et al., Proc. 26th ICRC, Salt Lake city, USA , HE2.2.09, 100 (1999) 

F. J. M Farley, and J. R. Storey, Proc. phys. Soc. , 996 (1954) 

H. Kojima et al., Proc. 27th ICRC, Hamburg, Germany , 3943 (2001) 

H. Kojima et al., Proc. 28th ICRC, Tsukuba, Japan , 3957 (2003) 

K. Nagashima et al., Nuovo Cimento Soc. Ital. Fis. , 695 (1989) 

K. Nagashima ct al., K. Fujimoto, and R.M. Jacklyn, J. Geophys. Res. , 17429 (1998) 

M. Aglietta et al., APJ , 501 (1996) 

T. Kifunc et al., J. Phys. G , 129 (1986) 

R. Clay et al., Proc. of 25th ICRC, Durban , 185 (1997) 

T. Antoni et al., ApJ , 687 (2004) 

D. Heck ot al.. Report FZKA6019 (Forschungszentrum Karlsruhe) , (1998) . 

tycho. usno . navy, mil / sidereal . lit ml 
R.M. Jacklyn, PASA , 425 (1986) 
H.R. Allan, ApJ. Lett. , 237 (1986) 
J. Linsley, Phys. Rev. Letts. , 1530 (1975) 

M. Aglietta et al, Proc. of 28th ICRC, Tsukuba, Japan , 183 (2003) 



7 



Amplitude(%) Pha5e(h) 
A/ =0.32 ±0.10 Tr = 21.3 ±1.0 
An = 0.56 ± 0.10 Tjj = 20.6 ± 0.7 



TABLE I: Parameters of the fit coefficients in equation(7) 
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FIG. 1: Schematic view of water Cherenkov detectors as a square array, and the electronic circuit. Particle tracks through lid 
(right hand- top), and tracks through walls( right hand-bottom) have also been shown. 
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FIG. 2: Distribution of events time-spacing. 
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FIG. 3: Event rates per half an hour as a function of atmospheric pressure. 
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FIG. 4: Distribution of mean 1/2-hr event rates before (squares) and after (points) the atmospheric correction for pressure. 
The line curve shows a fit by a Gaussian function. 




FIG. 5: Frequency of air showers vs. zenith angle. 
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FIG. 6: Observed sidereal time variation in intensity of the cosmic rays (points). The curve is the best fit to Eqn. (7) with the 
coefficients as listed in Table I. 



14 




FIG. 7: Celestial coordinate, C= direction of cosmic ray, M= direction of earth's motion, Z= zenith, P= direction of North 
pole. 
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FIG. 8: Distribution of air shower events vs. declination angle. 
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FIG. 9: The number of photons received by PMT in consequence of particle passing along different tracks entering through the 
lid( Top), and the walls( bottom). 



